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^ ABSTRACT 
O 

We calculate the theoretical expectation for the bulk motion of a large scale 
survey of the type recently carried out by Lauer and Postman. Included are 
the effects of survey geometry, errors in the distance measurements, clustering 
properties of the sample, and different assumed power spectra. We consider the 
power spectrum calculated from the IRAS-QDOT survey, as well as spectra from 
hot + cold and standard cold dark matter models. We find that sparse sampling 
and clustering can lead to an unexpectedly large bulk flow, even in a very deep 
survey. Our results suggest that the expected bulk motion is inconsistent with 
that reported by Lauer and Postman at the 90 — 94% confidence level. 



Subject headings: cosmology: theory - observation - galaxies: clustering - distances and 
redshifts 
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1. Introduction 

Recently, Lauer and Postman (1993) (LP) measured the bulk motion of a volume limited 
complete sample of all 119 Abell clusters with redshifts out to 15,000 km/s. They used a 
distance indicator based on brightest cluster galaxies that gives a typical distance error of 
16% for a single cluster. They report that the frame defined by their sample exhibited a bulk 
motion of 689 ±178 km/s with respect to the cosmic rest frame defined by the CMB. This 
value seems quite large given current expectations for the power spectrum on these scales. 
However, as we shall see below, the Abell clusters are a sparse sampling of the underlying 
peculiar velocity field, in the sense that the intracluster spacing is large compared to the 
scale on which the field has significant variations. This implies that the bulk motion of the 
Abell cluster sample is not necessarily indicative of the motion of the volume as a whole. 

In this letter, we study the theoretical expectations for a measurement of the type 
made by Lauer and Postman. Given a catalog of cluster positions and an assumed power 
spectrum, we calculate the covariance matrix for a measurement of bulk motion in the sample 
using a linear analysis based on the work of Kaiser (1989). We show that the covariance 
matrix is a sum of a noise term and the convolution of the velocity power spectrum with 
a window function. Given that the covariance matrix is roughly diagonal and isotropic, we 
can calculate the expectation value for the bulk flow as well as the probability for measuring 
a velocity greater than some reference value. We consider power spectra from the IRAS- 
QDOT survey (Feldman, Kaiser & Peacock 1993), the BBKS CDM model (Bardeen et al. 
1986) and from simulations of hot + cold dark matter (HCDM) (Klypin et al. 1993). For 
our analysis we use the Lauer and Postman sample as well as simulated cluster catalogs 
designed to explore the dependence of our results on the radius and clustering properties 
of the survey. It should be noted that we treat the power spectrum and cluster catalogs as 
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being independent. This allows us to avoid questions of whether the clustering properties of 
Abell clusters are representative of clusters as a whole (Efstathiou et al. 1992). In fact, we 
show that the high degree of clustering in the Abell cluster sample results in a significant 
increase of the expected bulk motion. Preliminary results of this analysis were presented 
elsewhere (Feldman & Watkins 1993). Our analysis can be contrasted with that of Strauss 
et al. (1993) who analyzed Monte-Carlo realizations of peculiar velocity data drawn from 
N-body simulations of specific theoretical models. 



2. Analysis 



The maximum likelihood solution for the uniform streaming motion Ui of a cluster 
sample (Dressier et al. 1987) is given by 

(7 -+- (7 

n J n 1 * 

where r n is the position of the nth cluster, S n is its estimated line-of-sight velocity, <r* is the 
ID velocity dispersion, o n is the estimated uncertainty in the line of sight peculiar velocity, 
and 

n u n ' * 

(T n is proportional to the distance and is denoted by o n = (/3/100)r n , where f3 is the per- 
centage error in the distance measure. For the LP survey f3 = 16%. In these equations and 
the equations to follow, repeated indices denote implicit sums. The estimated line-of-sight 
peculiar velocity S n is related to the true velocity by 

Sn = f'n,i v i( r n) ~\~ &n ? (3) 
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where e n is drawn from a Guassian with zero mean and variance a 2 e = a 2 n + a 2 . For the 
catalogs of interest Aij is nearly diagonal, with the off-diagonal terms being of order f 0% of 
the diagonal ones. 

Following Kaiser (1989) we write the uniform bulk velocity in terms of Eqs. (1), (2) and 

(3): 



Ui = U\* + u\ e) = A~ x £ r -f^v k (f n ) + A~ x J2 I tri ■ ( 4 ) 

n J n 1 * n n 1 * 

Since the two terms are statistically independent, the covariance matrix is 

R l3 =< U,U 3 >= + 4 £) • (5) 
The velocity term can be written as 

U\ V \r ) = f d 3 r W tl (f) Vj (f + r) , (6) 
where Wij is the tensor window function for the survey, given in fourier space by 

W ij (k) = A£ E%T% (7) 

n n ' * 

Then the velocity part of the covariance matrix is the convolution of a window function and 
the velocity power spectrum 

^ = j^f p v(k) , (8) 

where 

= Wu(k) w jm (k) h L . (9) 

The velocity power spectrum is 

P v (k)=< \ v (k)?>=^P(k) , (10) 
where P(k) is the density power spectrum. Like Wfj and thus are nearly diagonal. 
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As stated above, we consider power spectra from the IRAS-QDOT survey, the BBKS 
CDM (<r 8 = 1 fi/j = 0.5) model, and HCDM simulations (normalized to COBE quadrupole 
Q2 = 17//K). Since we anticipate our theoretical expectations to be low compared with the 
observed bulk flow, we will also include the IRAS-QDOT spectrum plus its la error in our 
study. We have corrected the IRAS-QDOT and HCDM spectra for the redshift distortion 
pointed out first by Kaiser (1987). 

The sparseness of the surveys we are considering will lead to features in the window 
function which will result in significant contributions to the bulk motion from scales much 
smaller than that of the survey. This is demonstrated in Fig. 1, where we show window 
functions for the LP and simulated catalogs in conjunction with the velocity power spectrum. 
Note that the window functions do not fall to zero outside of the central peak as they would 
if the volume were densely sampled. We also show a Gaussian exp( — k 2 R%), where is the 
effective depth of the survey. For volume limited surveys of this type, the effective depth of 
the survey is ~ 

Rmax/v^i where i? max is the volume radius. 
The noise term is 

4>^'E^)(^E^)H'. (") 

Thus the final covariance matrix [see eq. (5)] is 

Rii = J P v (%) + Ar> (12) 

Given that is nearly diagonal and isotropic, we can average over the three diagonal 
components and take 

a^^/S; a[ v) = y/^]3, (13) 

to be the ID variance for each of the three components of the bulk flow. Here a superscript 
v refers to the error free quantities. The probability for measuring a bulk flow of magnitude 
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V then becomes 

p(y)=Q 1/2 ^ex P (-y 2 /2- 1 2 ) • (14) 

The expectation value for the magnitude of the bulk flow in a given sample is 

A = r V P(V) dV = -^=a x ■ AW ee -4=<t[ v) . (15) 
Jo y / 2ir V27T 

Given one can also use P(V) to calculate the probability for observing a bulk flow greater 
than a given magnitude. 

We construct three types of simulated cluster catalogs. Two with correlation length 
r = 23 and 14 h~ x Mpc, using a method similar to that of Postman et al. (1989) and 
Plionis et al. (1992), and a third where clusters are distributed randomly, i.e. r = 0. 
As usual we parametrize the uncertainty in our knowledge of the Hubble constant by h = 
if/ (100 km s _1 Mpc _1 ). In order to mimic the LP catalog, we have eliminated all clusters with 
Galactic latitude \b\ < 15° and applied the esc extinction law P(b) = dex [0.13(1 — esc | &| )] 
(Strauss et al. 1993). For the 15, 000 km/s samples, we chose 120 clusters as in the Lauer and 
Postman catalog. To construct catalogs with different radii we assume a constant density of 
clusters and scale the number of clusters accordingly. 

3. Results 

In Tables 1-3 we present our results for both the real and simulated cluster catalogs 
for the four power spectra we consider. For comparison with our results, we are interested 
in the uncorrected measurement of the bulk flow Vlp. This value is obtained by taking 
the magnitude of the bulk flow vector reported by Lauer and Postman in the cosmic rest 
frame, namely Vlp = \F\ = |477 km/s, —142 km/s, 635 km/s| = 807 km/s. The bulk flow 
magnitude reported by LP, 689 ± 178 km/s, has been corrected for error bias. Note that this 
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'bias-free' magnitude is not the same as our 'error-free' magnitudes and they cannot 
be directly compared. 

In Table 1 we present the results of our analysis applied to the LP catalog. We show 
both the raw (A) and the error-free (A^) expectation values for the different power spectra. 
We also present the probability for measuring a bulk motion larger than that observed by 
LP, PiV > Vlp). Results for the IRAS-QDOT spectra should be taken as lower bounds, 
due to a the lack of reliable information for k < 0.02 (see Feldman et al. 1993). These 
results indicate that the power spectra we have considered are inconsistent with the Lauer 
and Postman observation at the 90 — 94% confidence level. 

In Tables 2 and 3 we present the results for the simulated cluster catalogs. The results 
shown here are averaged over 50 realizations. The statistical variations in A and A^ were 
found to be of order 5% and so are neglected. In Table 2 we present the results for both 
A and A^ for the three classes of clustering properties and measurement errors (see Figure 
2a-c). We see that the results of the highly clustered simulated catalogs are very similar 
to those from the LP catalog, suggesting that they have similar properties. Here we used 
f3 = 16%, as reported by LP. The clustered catalogs enhance the bulk flow relative to random 
catalogs by 10 — 20% (see Figure 2c). 

Lauer and Postman, in collaboration with Strauss (Strauss 1993 private communica- 
tion), are planning to undertake a project whereby they will survey all Abell clusters within 
24, 000 km/s with improved measurement errors of f3 = 10%. We have extended our analysis 
to include larger simulated surveys with the same cluster density. In Table 3 we show the 
results for a 25,000 km/s survey for f3 = 10% (also see Figure 2a-c). These results indicate 
that the bulk flow in this extended survey is expected to be ~ 250 km/s. 
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For comparison, in Tables 2 and 3 and Figure 2d we show the results for the 'true' 
expected bulk velocity of the frame in the limit of infinite number of randomly placed 
sample objects for the HCDM and CDM power spectra. 

4. Conclusions 

We have presented a formalism to calculate the theoretical expectation for bulk flow 
in large scale surveys as a function of the geometry of the survey, its clustering properties 
and the assumed power spectrum. We apply the formalism to the LP survey as well as 
to simulated cluster catalogs. We find that there are two effects that significantly enhance 
the bulk motion for a catalog, namely sparseness and clustering. Sparse sampling of the 
velocity field leads to contributions from the velocity power spectrum on smaller scales than 
the survey. For a LP type survey this can result in a large enhancement over the expected 
motion of the volume as a whole. Clustering of the sample leads to even more enhancement, 
since it effectively increases the sparseness. For the clustering in the Abell catalog this results 
in an additional fO — 20% enhancement. 

Our results indicate that the power spectra we have considered are inconsistent with 
the LP measurement of bulk flow at the 90 — 94% confidence level. While this level of 
inconsistency is provocative, it does not indicate, as yet, a major problem with current 
theories of large scale structure formation. Improvements in the accuracy of the distance 
indicator should help clarify this issue. 

We have also explored the convergence of the survey frame to that of the CMB frame 
as a function of radius using simulated cluster catalogs. Our results indicate that an Abell 
cluster survey of radius 25,000 km/s with a distance measure accurate to 10% will exhibit 
a bulk motion of order 250 km/s, although this is somewhat uncertain due to a lack of 
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information about the power spectrum on these scales. Extending the survey to this radius 
will provide interesting information about scales where the uncertainty in our knowledge of 
the power spectrum is large. 

Acknowlegements: We are grateful to Nick Kaiser, Gus Evrard, and Andrew Jaffe for 
valuable comments. We also thank Tod Lauer and Mark Postman, and Michael Strauss for 
answering our questions regarding their work. HAF was supported in part by the National 
Science Foundation grant NSF-PHY-92-96020. RW was supported in part by NASA grant 
NAGW-2802. 



- 10 - 
REFERENCES 

Bardeen, J. M., Bond, J. R., Kaiser, N. & Szalay, A. S. 1986, ApJ, 304, 15 

Dressier, A., Faber, S. M., Burstein, D. Davies, R. L., Lynden-Bell, D., Terlevich, R. J. & 
Wagner, G. 1987, ApJ, 313, 42 

Efstathiou, G. Dalton, G. B., Sutherland, W. J. & Maddox, S. J. 1992, MNRAS, 257, 125 

Feldman, H. A., Kaiser, N. & Peacock J. A. 1993, ApJ in press 

Feldman, H. A. & Watkins, R. 1993, proceedings of IAP Cosmic Velocity Fields, Ed. F. 
Bouchet & M. Lacheieze-Rey, in press 

Kaiser, N. 1987, MNRAS, 227, 1 

Kaiser, N. 1988, MNRAS, 231, 149 

Klypin, A., Holtzman, J., Primack, J. & Regos, E. 1993, ApJ, 413, P48 

Lauer, T. & Postman, N. 1993 preprint, LP 

Plionis, M., Valdarnini, R. & Yi-Peng, J. 1992, ApJ, 398, 12 

Postman, M., Spergel, D. N., Sutin, B., & Juszkiewicz, R. 1989, ApJ, 346, 588 

Strauss, M., Cen, R. & Ostriker, J. P. 1993, proceedings of IAP Cosmic Velocity Fields, Ed. 
F. Bouchet & M. Lacheieze-Rey, in press 



This preprint was prepared with the AAS IATj^X macros v3.0. 



Table 1: LP Catalog (V LP = 807 km/s) 
A and are given in km/s. 
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P(V > Vlp) 
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470 
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0.06 


IRAS + a 
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0.10 


CDM 


240 
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0.07 


MDM 
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0.07 



Table 2: i? max = 15000 km/s; fi = 16% 



A, AH and AW are given in km/s, r values are in h 1 Mpc. 
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Table 3: i? max = 25000 km/s; fi = 10% 



A, AH and AW are given in km/s, r values are in h 1 Mpc. 
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Figure 1 



Fig. 1. — Fig. 1) In [a] we present the redshift corrected power spectra used in the analysis. 
In [b] we see the normalized tensor window for the highly clustered and random catalogs, 
the one for the LP catalogue and a guassian window function exp( — k 2 R%). 
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Figure 2 

Fig. 2. — Fig. 2) In [a] we show the expectation values as a function of survey radius 
for f3 = 10% (A) and 0% (error- free A^) for the IRAS-QDOT power spectrum. In [b] 
we show the error-free expectation values (A^) for the four power spectra we used. In 
[c] we show the expectation values for {3 = 10% (A) and the error-free (A^) for random, 
moderately clustered and highly clustered (labeled 0, 14, 23 respectively) simulations for the 
IRAS-QDOT power spectrum. In [d] we show the results for the limit of an infinite number 
of clusters for the CDM and HCDM power spectra (A^). 



